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ABSTRACT 

Aims. We analyzed the regions of star formation in the spiral galaxies NGC 3184 and NGC 3938 from archive images 
for a wide range of wavelengths (NUV from GALEX, Ha from JKT and KPNO, 8 and 24 /im from Spitzer, and CO 
from BIMA). 

Methods. We used the Clump Find Algorithm to extract the properties of the star-forming tracers, identifiable as 
emission regions at each wavelength. 

Results. We obtained a power-law relation between the luminosity and the emission region volume that scales as 
expected, L tx V, for the Ha and NUV emission, but the luminosity varies far more rapidly with the volume for the 
dust (8 and 24 /jm) and molecular gas emitting regions in CO. This is interpreted as a change in the emissivity with 
the size of the cloud, either by an augmentation of the overall density or due to the presence of high density clumps, 
with high local emissivity coefficients. Although the clumpy nature of molecular gas may be unsurprising, the clumpy 
nature of mid-infrared emission regions, which could be explained by newly formed high to intermediate mass stars 
embedded in the dust providing the heating, is clearly revealed in both galaxies. 
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1. Introduction 

Surveys completed with both the Spitzer Space Telescope 
(Werner et al. 2004) and the Galaxy Evolution Explorer 
(GALEX; Martin et al. 2005) provided morphologically de- 
tailed IR and UV images of galaxies. Likewise, surveys 
such as the Berkeley Illinois Maryland Association (BIMA) 
Survey of Nearby Galaxies (SONG; Heifer et al. 2003) pro- 
vided us with aperture synthesis CO mapping surveys of 
galaxies. These surveys are enabling multiwavelength stud- 
ies of star-formation activity with adequate resolution and 
sensitivity (e.g. Calzetti et al. 2005, 2007; Kennicutt et al. 
2007; Thilker et al. 2007), and the cross calibration of star- 
formation rate (SFR) indicators at different wavelengths 
(see Calzetti et al. 2007). 

Since the work by Kennicutt (1988), studies of HII re- 
gions in spiral galaxies by means of their Ha emission have 
demonstrated that the logarithm of the diameter of the 
H II region and the logarithm of the Ha luminosity scale 
roughly with a value of one third, i.e. linearly with the 
volume, which is the slope expected for constant density 
radiation-bounded nebulae. This trend is found in studies 
of H II regions of several spiral galaxies, such as NGC 7331 
(Marcelin et al. 1994), NGC 3992 (Cepa & Beckman 1989), 
and NGC 4321 (Cepa & Beckman 1990), as well as in irreg- 
ular galaxies such as NGC 6822 (Hodge et al. 1989), and 
IC 10 (Hodge & Lee 1990). 

We present a study of star-formation regions in the 
spiral face-on galaxies NGC 3184 and NGC 3938 using im- 
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ages at various wavelengths (NUV from GALEX, Ha from 
JKT and KPNO, 8 and 24 /mi from Spitzer, and CO from 
BIMA). Our aim is to explore the possible slope values of 
the relation L v oc V$ for the collection of individual gaseous 
or dusty star-forming regions of volume V and luminosity 
L v detected at each frequency bandwidth v. 

This work is a phenomenological study to provide clues 
in understanding the different processes affecting the ob- 
served beta slope. 



2. Multiwavelength Images of NGC 3184 and NGC 
3938 

The galaxy images were taken from archive data. We 
used the Spitzer Infrared Nearby Galaxy Survey (SINGS; 
Kennicutt et al. 2003) database to select candidate spi- 
ral galaxies oriented nearly face-on, which are the most 
suitable for our study because optical and UV data are 
less affected by extinction. From a short list of candidates, 
we selected NGC 3184 and NGC 3938, both with avail- 
able 8 /im Infrared Array Camera (IRAC; Fazio et al. 
2004) images and 24 /im Multiband Imaging Photometer 
(MIPS; Rieke et al. 2004) images from Spitzer with a res- 
olution of 0.75 " pixel -1 . For both galaxies, CO maps at 
2.6 mm (J:l— >0) are also available with a resolution of 
0.997" pixcP 1 x 1.009 "pixel" 1 from BIMA SONG (Heifer 
et al. 2003). The Ha image of NGC 3184 with a resolu- 
tion of 0.234" pixel -1 is from the archive of the Jakobus 
Kapteyn Telescope of the Isaac Newton telescopes (Knapen 
et al. 2004), and that of NGC 3938 was obtained from the 
NASA/IPAC Extragalactic Database (NED) taken at the 
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4. L u OC 




Fig. 1. NGC 3184 . Comparison of the original images (left) and 
those obtained with the Clumpfind2D algorithm (right) at Ha 
(top) and 8 /xm (bottom). Both images cover an 8' x 8 'region, 
while in the bottom right panel the solid line square is the area 
studied in CO (5.7' x 5.6') and the dashed square is the 24 /jm 
image (4.5 ' x 4.5 '). 



2.1m telescope at Kitt Peak National Observatory (KPNO) 
of NO AO (National Optical Astronomy Observatory) with 
a resolution of 0.303 " pixel -1 x 0.302 " pixel" 1 . The 2271 A 
(NUV) image was obtained by GALEX (Martin et al. 2005), 
and has a resolution of 1.513" pixel -1 x 1.514" pixel -1 , and 
covers a bandwidth of 1750-2800 A. 



3. Clump Finding Analysis 

To be able to study the clumpy structure of the galaxies 
at each wavelength, we used the Clump Finding algorithm 
developed by Williams et al. (1994) to study structure in 
molecular clouds. We have used a 2D version of the rou- 
tine, Clumpfind2D, originally in Fortran and then adapted 
to IDL. This algorithm defines different clumps using the 
"friends-to-friends" routine, which groups pixels with some 
sort of connectivity among them, and inside a specific range 
of intensity. The final result of the Clumpfind2D routine is 
a list of the different clumps, their position, and their area 
in pixels. 

Galaxy images were preprocessed. Foreground stars 
were removed from the original galaxy images using stan- 
dard IRAF routines. The areas explored for clumps in each 
galaxy correspond to the minimum area in which the galaxy 
is observed at all wavelengths. For NGC 3184, the areas 
studied are Ha 8 ' x 8 ', 8 /im 8 ' x 8 ', 24 /mi 4.5 ' x 4.5 ', 
and CO 5.7' x 5.6'(c.f. bottom right panel in Fig. 1). For 
NGC 3938 , a similar area of 5.7' x 5.6' was studied at all 
wavelengths, NUV, Ha, 8 and 24 /xm, and CO. 

As an example of how the morphology is well repro- 
duced, in Fig. 1 we show the original image of NGC 3184 
at Ha (top) and 8 /jm (bottom) with the corresponding 
images obtained by Clumpfind2D (right panel). 



The flux-calibrated original image and the clump or 
emission-region parameters obtained with ClumpFind2D 
allow us to calculate the luminosity of each region emit- 
ted at a particular wavelength. To calculate emission-region 
luminosities, we assumed a distance of 11.5±2.0 Mpc for 
NGC 3184 and of 11.6±0.8 Mpc for NGC 3938 , taken from 
NED. Different corrections were applied to the calibrated 
fluxes depending on the wavelength, for NUV, a Milky Way 
(Cardelli et al. 1989) and internal extinction (Buat et al. 

2005) corrections; for Ha, a Milky Way extinction (Cardelli 
et al. 1989) and contribution of N [n] (Calzetti et al. 2007) 
corrections; and for 8 /xm, a Milky Way extinction correc- 
tion (Indcbctouw et al. 2005). No corrections were made 
for CO and 24 /xm fluxes. Luminosities were calculated at 
the corresponding frequencies of Ha, 8 and 24 /xm, CO 2.6 
mm, and NUV (2271 A) if available. We therefore calculate 
luminosities of regions of emitting gas (ionized or molecu- 
lar) or dust at 8 and 24 /xm, which have been identified as 
characteristic tracers of star-forming regions. 

We noticed a relation between the areas of the emis- 
sion regions and their luminosities for a given wavelength. 
Assuming a spherical geometry, the volume was found 
based on the equivalent radius obtained from the original 
emission-region area. The relations between the multiwave- 
length luminosities L v of the clumps and their correspond- 
ing volumes V (defined to be 47ri? 3 /3, where R is the mean 
radius of the emission region) obtained for NGC 3184 and 
NGC 3938 are of the form L v oc V 13 . 

To visualize the difference between all of star-formation 
tracers, we present them in different panels on Fig. 2 and 
plotted together on Fig. 3. In both galaxies, the disper- 
sion in Ha is comparable to previous results (e.g. Kennicutt 
1988). For CO and NUV, the dispersion is comparable with 
that in Ha, while for both 8 and 24 /xm the dispersion is 
smaller. 

For NGC 3938, the slope values (/?) for both Ha and 
NUV are close to 1, while for the other tracers the slopes 
are much steeper (see the right panel of Fig. 2). On the 
other hand, for NGC 3184 the slope for Ha does not become 
close to 1, but the other luminosities follow the same trend 
with volume as that found in NGC 3938 . For both galaxies, 
the dust tracers render the most similar (3 values, being 
identical for 24 /im. The steepest value is obtained for the 
two galaxies with the CO clumps. 

5. Discussion 

5.1. Star-formation tracers 

The various wavelengths are related to the star-formation 
process in different ways. The NUV emission originates in 
regions populated by massive, young stars (ages < 10 8 yr). 
This radiation is affected significantly by extinction, but af- 
ter the necessary corrections have been made, its intensity 
relates directly to the actual SFR (Iglesias-Paramo et al. 

2006) . The Ha emission is the principal recombination line 
used to trace the transparent star-forming HII regions, 
where the presence of dust is negligible and molecular gas 
has been wiped out (Osterbrock 1989). There has been a 
discussion about the problems of using the 8 /xm emission 
dominated by PAHs as a tracer of star formation (SF), and 
no final conclusion has been reached (Calzetti et al. 2007; 
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Fig. 2. Relation between the luminosity and volume of the emission region obtained for different star formation tracers in NGC 
3184 and NGC 3938. 



Wu et al. 2005). Since PAHs have been found to be associ- 
ated with very young stars, it is likely that its emission is a 
reliable SF tracer. The relation between SF and emission at 
24 fxm due to a hot dust phase has been investigated sev- 
eral times (see Helou et al. 2000), and is found to originate 
in very small dust grains (VSG) in regions surrounding the 
ionized gas. Finally, CO traces the molecular clouds inside 
which the process of SF occurs (molecular cores) or will 
occur (Kennicutt et al. 2007). 

5.2. On the value of the slope (3 

The luminosity of an emitting cloud of gas (or even stars) 
is to a first approximation given by L v oc J v e u dV ', where 
the integral is computed over the entire volume of the cloud 
V and e„ is the local emissivity at frequency v. If all pa- 
rameters that control e„ are constant within the volume, 
the expression reduces simply to L v oc V. If the shape of 
the cloud is globally spherical, then L u oc V oc R 3 . In this 
case, we expect a slope of /? = 1, which is observed in the 
Ha image of NGC 3938 and in the NUV image of the same 
galaxy. 

5.2.1. The case of (3 > 1 

To derive a value of (3 that differs from one, one of the 
previous hypotheses must be changed: 



— The shape must not be spherical: the clouds could be 
instead elongated along the line of sight, and this elonga- 
tion should increase with an increasing size of the cloud. 
This hypothesis, however, appears very unrealistic. 

— The local emissivity e v could vary with the size of the 
cloud. This is possible if any parameter to which the 
emissivity is related varies as well, for example the den- 
sity of the emitting gas. Most of the emissivity coeffi- 
cients are related to the square of the gas density. If 
this is the case, we could obtain a slope larger than 
one if the mean density of the gas increases with the 
volume of the cloud. This could be due to a global in- 
crease in the density, or to the presence of high density 
clumps, with high local emissivity coefficients. For these 
clumpy regions, the resulting luminosity can be decom- 
posed into two contributions: L v oc e°(Vr — V c ) + e£V c , 
where Vr and V c are the total volume of the cloud and 
the volume of the high density clumps, respectively. 
The emissivity coefficient of the background low den- 
sity gas is e° and the clump emissivity coefficient is 
t c v . If the volume of the clumps is related to the to- 
tal volume of the cloud according to V c — kVr, where 
k is similar to a filling factor, we can describe the lu- 
minosity by L oc [e°(l — k) + e^k]Vr- To observe any 
effect of the clumps, the contribution of both phases 
of the gas should have the same order of magnitude: 
e" (1 — k) ~ e c v k. Since the clumps have higher emissivity 
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coefficients (e£ > e°) and a smaller relative volume (A; 
small), one of these two parameters must increase with 
the size of the cloud to obtain a luminosity higher than 
that expected from pure volumetric effects (/? > 1). The 
first option is that the emissivity of the dense clumps 
should increase: we can always consider denser clumps 
as the total size of the cloud increases. The second op- 
tion is that the relative volume of the clumps increases, 
i.e. there are more or larger clumps in larger clouds. 
This would also explain the case of (3 < 1 for Ha in 
NGC 3184 (see below). 

— It is known that due to density inhomogencitics, molec- 
ular gas tends to form clumps. However, clumps found 
at mid-infrared (MIR) emission regions (which could be 
explained by newly formed high to intermediate mass 
stars embedded into the dust that heat and re-emit 
with a temperature dependence T 4 ) are probably due 
to a global increase in the density, or to the presence 
of high density clumps, with high local emissivity co- 
efficients. The clumpy nature of the dusty regions is 
clearly revealed in both galaxies. As shown in Fig. 1, 
MIR clumps closely follow the spiral arms and reach far 
larger volumes than the other tracers (see Figs. 2 and 
3). The ionizing regions detected by their Ha and NUV 
emission appear to span a smaller range in volume than 
dusty and molecular gas regions, which reach far larger 
volumes. 

5.2.2. The case of (3 < 1 

The same study of the hypothesis corresponding to (3 = 1 
must be completed to understand the possible situations 
leading to (3 < 1 (which was obtained for the Ha image of 
NGC 3184): 

— The shape of the cloud is not spherical: it must smoothly 
change from being spherical to oblate ellipsoidal, with 
a decreasing size in the direction of the line of sight. 
In the present case of face-on galaxies, this direction is 
perpendicular to the plane of the galaxy. The value of 
(3 would change from 1 to 2/3, as the luminosity would 
become proportional to the area. This is not what is 
actually observed, since the slope is rather constant for 
a large range of cloud sizes. 

— Another way of changing the shape of the emitting cloud 
is to take into account the optical depth along the line 
of sight, which is related to the cloud's size in this direc- 
tion. Increasing the size of the cloud would lead to high 
values of the optical depth, such that the most distant 
part of the cloud could not be seen. In this case, it would 
be the shape of the visible cloud that would change, 
from a naturally spherical shape to a more disk-shaped 
when observed face on. We would observe a type of sat- 
uration effect as the optical depth started to become 
important, which would again change the slope from 1 
to 2/3, something that is not observed. 

— The volume of the emitting cloud can also be reduced 
by removing its inner part. This may be the case if there 
were some places inside the cloud from which the gas 
does not emit. This could be explained if some clumps 
of gas are in a phase in which they do not emit. The 
presence of neutral or molecular high density insertions 
(high density clumps inserted on less dense surround- 
ings), which not only do not emit in Ha, but also con- 



ceal the gas behind them, can lead to a decrease of the 
emitting volume as the size of the cloud increases. Such 
insertions would also explain the (3 > 1 slope in the case 
of dust or CO maps as is observed (see Sect. 5.2.1). 
— We did not explore the possibility that emissivity de- 
creases with the size of the cloud, since we cannot iden- 
tify any physical process that could be responsible for 
this behavior. 

6. Conclusions 

We have studied star-forming regions in the face-on spiral 
galaxies NGC 3184 and NGC 3938 , which have emission 
detectable by a variety of different tracers: NUV, Ha, CO, 
and 8 and 24 /zm. By analyzing the multiwavelength lumi- 
nosities L v of star-forming emission regions and their corre- 
sponding volumes V, we have found a relation of the form 
L v oc V 13 . If we suppose that the luminosity is directly pro- 
portional to the volume, we should expect a value of (3 =1, 
which is indeed observed for the images in Ha and NUV in 
NGC 3938 . 

However, from the CO and 8 and 24 fim images, a value 
of P > 1 is derived for both NGC 3184 (j3 C o ~l-8) and 
NGC 3938 {(3 C o ~1.5). We conclude that this could be due 
to a change in the local emissivity e„ with the size of the 
cloud, either by an augmentation of the overall density or 
because of the presence of high density clumps, with high 
local emissivity coefficients. For both galaxies we measure 
similar values of (3 for the dust tracers, i.e. (3 ~1.2-1.3 for 8 
fim and 1.51 for 24 //m. From previous studies, the clumpy 
structure found for molecular gas and Ha was expected. 
The mid-infrared emission of dust (PAH and VSG) heated 
by newly formed high to intermediate mass stars is revealed 
as clumps as well throughout the galaxies. 

For the Ha image of NGC 3184, a slope of j3 <1 is 
obtained, which could be interpreted by the presence of 
neutral or molecular high density clumps, which not only 
do not emit in Ha but also conceal the gas behind them. 
This can lead to a decrease in the emitting volume as the 
size of the cloud increases. 

For NGC 3184 and NGC 3938 , the volumes found for 
ionizing regions (Ha and NUV emission) appear to span a 
smaller range than that found for dusty (8 and 24 /zm) and 
molecular gas regions. 
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